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Abstract 

Renormalized vacuum expectation values of electromagnetic stress-energy ten- 
sor are calculated in the background spherically-symmetrical metric of the worm- 
hole's topology. Covariant geodesic point separation method of regularization is 
used. Violation of the weak energy condition at the throat of wormhole takes 
place for geometry sufficiently close to that of infinitely long wormhole of constant 
radius irrespectively of the detailed form of metric. This is an argument in favour 
of possibility of existence of selfconsistent wormhole in empty space maintained by 
vacuum field fiuctuations in the wormhole's background. 
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1. Introduction. In 1988 Kip Tliorne witli co-workers have studied properties 

of static spherically-symmetrical traversable wormhole. To exist, such the wormhole 
should be threaded by material with rather unusual properties. In particular, radial 
pressure of the material should exceed it's density both locally at the throat [jl| and 
integrally along radial direction [^. That is, weak energy condition (WEC) 0] at the 
throat and averaged weak energy condition (AWEC) should be violated. Where can 
we find such the material? In Casimir vacuum between conducting spherical plates 
surrounding the throat was considered as a kind of such the material since it possesses 
required properties. In the given note another possibility is studied: self-maintained 
wormhole threaded by electromagnetic field vacuum. Covariantly renormalized vacuum 
expectation values (VEV's) of electromagnetic stress-energy tensor are found to violate 
WEC at the throat of the wormhole if it's typical length (proper radial distance at which 
variation of radius becomes comparable to the radius itself) is large as compared with 
it's radius. As for AWEC, check of it requires more complicated calculation and will be 
given in forthcoming paper. 

2. Calculation. In calculation it is convenient to use instead of radius r the radial 
distance p, so interval will take the form 

ds"^ = r'^{p)dn'^ + dp^ - exp (2$)dt^ dfi^ = dO'^ + sin^^ (1) 

In the case of wormhole $ is everywhere finite Consider electromagnetic field in this 
metric background, F^^ = d^A^, — d^A^, in the gauge At = 0. There exist transversal 
electric TE and transversal magnetic TM modes with components of vector potential of 
the form 



TE ■.{A^,Ag,A^) = (0,^, -Ygstne)R (2) 

sm 

TM : {Ap,Ag,A^) = (1(1 + l)^Y, R,Ye, R^^^^ (3) 

Here Y = Y^OA), R = Rmip), oo = u^u Ye = dY/dO, Y^ = dY/d<j), Rp = dR/dp. 
Radial functions obey eigenvalue equation 

- (^expi^)-^J R + l{l + lf-^^^^R = u;'R (4) 
with boundary conditions on conducting spherical plates T (if any) 

TE: R\r = 0; TM : Rp\r = 0. (5) 

It is apparently conformally invariant. 

Substitute appropriately normalized solutions for it into split-regularized form of 
stress-energy tensor 
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x—^x 4 / 

(6) 

and sum over eigenmodes (over n, I, m and polarization a = TE or TM). Here notation 
F^v{x)x means components of tensor F^y{x) transported in parallel way from point x to 
X along the geodesic connecting both points. It is convenient to split the point only in 
coordinate p; then parallel transport defined by vanishing the covariant derivative means 
simply constancy of "physical" components, 



^ui>...A == \gtitiguu---g\\\ ^'"^A^u... lambda = COnst, (7) 



def 

of any transported in such the way tensor A. Notations coincide with those of ref.[|l[ 
and are convenient in what follows. Normalization of radial functions should enable 
contribution of each eigenmode to energy 



- J ng^'^dS, = -j T\g^/Hpded(^ [g ^ -det||^?^,||) (8) 

to take on vacuum value %uj/2 (this is just covariant component of vector, as it should 
be for energy quantum mechanically conjugated to contravariant time). By properties 
of spherical functions summation over m eliminates angle dependence. The resulting 
vacuum expectation values of components of stress-energy tensor will be denoted just 
as components themselves; this will not lead to any confusion. So we have for nonzero 
components 

T- = lim y ^ + 1/2 RpRp + uj'RRexp (-$ - $0 + RR'-^ 

" ™ n,i,h™ Stt^'^^ / exp {-<!>) dp 

T = y V ^ + 1/2 RpR, + iu'RRexp (-$ - $0 - RR^-^ 

n,ikTM J R^ exp (-$)rfp ^ > 

^ 1 + 1/2 RR}^ 
Tee = S E ^ 



n,Z,TE,TM 



Snurr J R"^ exp {—^)dp 



Here e = \p — p\, R = R{p), r = r{p) and hereafter the components of any tensor 
obtained by interchange 9 ^ cf) are the same and will not be written out. 

Due to formal conformal invariance of electromagnetic field it is natural to temporar- 
ily change variable p ^ z via dz = exp {—^)dp and to consider the sums 



E + o) E - (= 5'radial) 

I ^ n ^ 

I ^ n 
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1 p p 

E(^ + o)^(^ + i)E— = ^a(=W) (10) 

I ^ n ^ 

I ^ n ^ 

Then = f/^i and, by eq.(D , S{ = -U^2 + r'^ exp (2$)S'a so that 



exp (-$) / exp^ exp($)\ exp (-$ - $) 

OTTifi = 6a ^ \ z H Z ty^(5-^) 

^ exp (-$) exp ($) exp('l) ^^ ^ exp(-$-'l) ^^ 

87ripp = i)a 5-: z H Z Uz{z~z) (11) 

SttT^-^ = Sa 



Computation of sums over n in the above expressions reduces to that of a Green 
function by the following trick. We introduce 

^EorM^^EorM^^^^^^) _ ^ G = + (12) 

n,TE or TM ^ + ^ 

through which U and S'a are expressible: 

r°°2/ + l , _^2/ + l 



u = Y.r—Gdt^t—Gdt 

= ^^^l{l + l)Gdt (13) 

In turn, G^, G^ are some Green functions: 



- Gfi^ + l^t' + l{l + 1)^^^^^) G^-'^ = 5{z - z') 
G^\z(^T = 0, G^lzgr = 
These eqs. are then solved perturbatively. Put 



(14) 



exp (2$) exp (2$ 



oj 



r 
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V, (15) 

'0 

where index denote values of considered quantities at the throat (at p = 0), and expand 
over V. 

We are interested here in empty space and thus shift F to infinity; then it proves that 
G^ = G^ = G, G = 2G. In zero approximation {V = 0) 
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G(o) = ^exp(-m|^-z|), ^ + l{l + 1)^^^^^. (16) 
2m Tq 

In the first (linear in V) approximation 

G = 2&'\z,S)-2l{l + l) J &'\z,y)V{y)G^'\y,z)dy (17) 

Calculating so, we are faced with integrosums over /, dt. There are two basic ones, and 
these arise already in zero approximation. In this approximation, substituting variables 
t = exp (^Q)rQ^[l{l + 1)]^/^ sinh (y9 and then ecoshy) = r^q we get with the help of (0), 
(ll) and (0) 

-/i(e/ro) 

2vr\^ T^P I = I -/i(e/ro)-/2(e/ro) | (1^ 



where 




and function h{q) is regular at g = 0. Carefully expanding it under the integral sign in 
Taylor series we get 



1/^(0) 7r/^'(0) 1 



h"{q) - h"{0)e{M - q) 



TT h"\<d) _ l_ 
48 



°° /i"(g) - h"{Q)e{M - q) 



hit) 



12 e 48 [Jo q 
2 /l(0) TT /i'(0) 1 /l"(0) TT /i'"(0) 



(^In 


2M 


+ 




e 


2M 


31\ 






e 





3 e4 
1 

+ 24 



4 e3 



2 e2 



+ 



12 e 



For /;,(g) we have 



(20) 



h{q) 



00 , Is 

4 /(?) 



00 



rr(g)-4g/'(g) + 6/(g). 



1 



(21) 



The /(g) is also regular at g = 0. Thus /i(0) = 6/(0), /i'(0) = 2/'(0), h"{q) = gV"'(9), 
/i'"(0) = 0, h"{0) = 2/"(0). In /(g) we expand exponential denoting / + i = 7: 
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exp 
F = 



-g 7 



7 




exp (-g7) 



1 1 

^ 12873 



+ 0(7' 



(22) 



This is enough to express required derivatives of /(g) at zero in terms of those of el- 
ementary sum g{q) = gX]exp(— 57). Lengthy also elementary calculation ultimately 



gives 











Tpp 


\'-\ 






rcg \ 



4 2 



^,4 



(23) 



3 e2 60 e 
3e^+ ^(1^7-1) 
-^(Inf-I)/ 

for regularized stress-energy tensor in zero approximation (for space-time taken as direct 
product of sphere and Minkowsky plain). The nonlogarithmic 0{h^^) contribution in 
Ji, I2 cannot be accurately estimated within our elementary approach and is simply 
included into logarithm thus resulting in the appearance of effective cut off L = kr^ with 
some numerical value k of the order of unity. 

Now one renormalizes gravity action by adding to it (and subtracting from matter 
action) some functional Wdiv so that 



rpflU 

-'div 
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-1/2 



div 



S9 



(24) 



cancels divergences in matter T'*'' This results in renormalization of cosmological 
constant, Einstein gravity constant and Weyl tensor invariant term in the action, 
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(25) 



where C^^x^ is conformal Weyl tensor. The full expression for was given by Chris- 
tensen . Nontrivial here is appearance of finite part of leading to nonzero trace of 
renormalized stress-energy tensor. 



ren ree div ' 



(26) 



so that 



rpfi 

fj,,ren 



[-ISR^^Xu^R'"'^'^ + SSR^^R'"' - 25i?2 -18nR 



28807r2 V — ■ ^^^^M— ' — ---V (27) 

Christensen's result depends locally on Riemannian and metric tensors and on e'^, the 
vector tangential at x to geodesic connecting points x and x', with length e, the length of 
geodesic. In our case e'^ = (0, 0, e, 0). In zero approximation (space-time=spherex plain) 
the Riemannian tensor has the only (up to index permutations) nonzero component 
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Rnigi = To • Substitute these input values into Christensen's formula and subtract the 



resulting T-^^ from our answer 



for Ti^". 



All the divergent terms, e ^, e ^ and Ine 



ones thus get cancelled. On the contrary, their cancellation is a useful check of our 
calculation. In logarithms the UV regulator e is substituted by infrared one A. The 
finite part of T^^ turns out to be nonzero in this simple case of geometry only for one, 
pp-component: 



2nX 



(0) 



pp 



(28) 



div, finite \ 

i.e. it is completely consists of anomaly. The resulting renormalized tensor takes the 
form 











Tpp 


1"' = ( 






ren \ 



^In^ 

60 1. 
60 J. 

^ In + ^ 

60 L ^ 120 



(29) 



The p — t symmetry inherent in this simple case is restored just due to anomaly. 

We now consider deviation of $, r from constants. Since we cannot take into account 
such deviation precisely, we need some parameter to expand over it. This parameter may 
be the ratio of wormhole's radius r(0) to the typical length - the proper radial distance 
p at which the variation r(p) — r(0) is comparable to r(0). In practice, it is equivalent 
to expansion over the full number of differentiations ^. For example, $" and are 
considered in such the expansion on equal footing. At the throat odd derivatives vanish 
due to symmetry, and the nearest correction to zero approximation will be given by 
r"-terms arising in linear approximation over deviation of $, r from constants. This will 
lead to the following three types of contributions to T^^^ at the throat (p = 0). 

First, those due to dependence on e of metric factors in (0) taken at x; these follow 
by expanding in Taylor series in e. At the throat we have: 



exp ($) 



exp ($0 
1 



l + 2'^oe' + 0(e^) 



r-o+2^o'e^ + 0(e^) 



(30) 



(here = r"(0) etc.). 

Second, those connected with changing expression of |z — 2| = | /o exp (— $)(ip| enter- 
ing already found formulas for T^'^g ^^'^ in terms of geodesic length e; to take this changing 
into account one should substitute in these formulas e by 



l-^<^>oe^ + 0(e^) 
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Third, contribution coming from linear in V correction to Green function G. Here 
we simply substitute Taylor series for V around y = into perturbative expansion (|1^ 
for G and integrate over y; due to exponents integration for each Taylor term is finite. 
Appearing integrosums reduce to that of above considered type I2, but now q^'^h is less 
singular at g = and/or it gets multiplied by positive power of e. 

These contributions are collected in the resulting total correction TI:Z^^\ respectively: 



pp 



(1) 



reg 



+ 



2f 

4^ + 2f 
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3 £2 
8 r'o/ro 
3 







18 rg ^ 



2^1 

3 e2 

3 £2 

sio _ ±'[0. I j_' 

3 e2 36 ~^ 36 



18 I 



' 3 e2 973- 
. (J)" 1 (J)" 

^2 



9 r2 



V 



4 1-;,' 



(31) 



With the same accuracy we get for components of Riemannian tensor at the throat 
(up to index permutations and 9 ^ (p): 



R 



_o 



, R, 



ipip 



(32) 



Others are zero in this approximation. Using Christensen's formula for T^^j^ with this 
input and subtracting it from Tj^^^ cancels all the singularities, as it should; besides that, 
we have for correction to finite part: 



(1) 



/- 



1 n 



pp 



div, finite 




73" 



V 



12 r-o 

1 r" 
— 9 
12 71 



72 r2 



(33) 



With taking into account all the contributions we ultimately get for total renormalized 
tensor in the considered approximation: 



2n' 



pp 



(o)+(i) 



/ 



^In^ 
1 In^ 







60r, 
1 

60^ 



lnA + 



_5_^ 

36 r, 



36 ~^ 72 / 



'0 



(34) 



3. Discussion. Taking the difference between radial pressure r = —Tpp and energy 
density g = obtained we find it to be positive at the throat: 



I n 
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This is quite nontrivial result, since generally we might get arbitrary combination of 
$Q, Tq/tq, not necessarily simply Tq/tq with positive coefficient. Thus, WEC is vio- 
lated at the throat, at least on the level of second derivatives: not only LHS of gravity 
eqs. (Einstein tensor) implies r > g, but also this condition is fulffiled by RHS (elec- 
tromagnetic stress-energy tensor) for given topology. As for the derivatives of larger 
order, their contribution to r — ^ is not a'priori fixed in sign, but the fact that required 
sign is achieved already on the level of second derivatives, means that contribution of 
higher derivatives required for further adjustment of given value need not be larger then 
that of ffist ones, which is normal condition of applicability of perturbation theory and 
regularity of possible solution. 

Large logarithms do not contribute r — ^ in the given approximation. This is because 
corrections to the structure 

9~'^'t^ I C,^,^C^^'-g"'d'x (36) 

entering T^^ as coefficient at logarithm, prove to occur only at the level of four differen- 
tiations. To see effect of the fourth derivatives we linearly expand (|36D . Einstein eqs. at 
the throat read: 



Rt 



R 



VP 



\R9it 



1 



R9, 



pp 



Rao 



\R9e 



1 
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27r2 
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1 
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27r2 \ 


V 604 


< + 


'_0_ 

ro 






. 1 


I N 


1+r 


- 


27r2 1 
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1 + 24 



5 rl 
364 



_2_ 

^0 



+ 



1204 



14 
36 4 



+ 



72 rl , 



(37) 



Here iV 



In^ 



SttG, G is Newtonian gravity constant. Not shown are bilinear in 
second derivatives terms multiplying which in precise calculations should be taken 
into account on equal footing with fourth derivatives. From pp - equation the estimate 
for wormhole radius follows, 



/240vr2A2\ 

(we have taken into account that itself depends on tq through L 
it follows from two other eqs. that 



krQ, k 



r" 



„2 ' 



' 



(38) 
Then 

(39) 



In other words, the length of the wormhole is of the order of it's radius. Thus, approxi- 
mation of long wormhole seems to be not relevant to selfconsistent solution, but we can 
hope to approach qualitatively physical wormhole starting from the sufficiently long one. 
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As for the range of possible values of ro, it follows from the experimental limitations 
on A. The latter arise in connection with that the induced C^-term should violate 
experimentally observed planetary motion, especially that of mercury The C^-term 
enters action with coefficient just of the order of Tq/k (the Tq itself as given by (^) is 
defined by this coefficient), i.e. with coefficient ~ relative to Einstein i?-term. So the 
effect of C^-term has relative value like exp (— rinerc/''^o) 0- Here rmerc is typical radius 
of mercury orbit. So the requirement for this correction not to exceed experimental 
accuracy leads to tq being at least 1.5 orders of magnitude smaller than the radius of 
mercury orbit, i.e. of the order of the radius of sun. As for the lower limit for tq, it can 
be estimated by taking A ~ 10^^ in Plank units, which corresponds to experimentally 
observable part of Universe. 

Thus, the range of experimentally acceptable (in indirect way) values of tq is quite 
large: it can be by half an order smaller than the Plank length (though in this region 
the concept of classical gravity may be wrong) and it can be as large as radius of sun. 

Finally, these estimates do not take into account the influence of another flelds, of 
which the most essential may be that of massless particles (e.g. the neutrino one) due 
to occurrence of large logarithms. Corresponding work is in progress now. 

I am grateful to prof. A. Niemi, S. Yngve and personnel of Institute of Theoretical 
Physics at Uppsala University for warm hospitality and support during the work on this 
paper. 
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